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Abstract 

We identify quark and gluon helicity-fiip distributions defined between nu- 
cleon states of unequal momenta. The evolution of these distributions with 
change of renormalization scale is calculated in the leading-logarithmic ap- 
proximation. The helicity-flip gluon distributions do not mix with any quark 
distribution and are thus a unique signature of gluons in the nucleon. Their 
contribution to the generalized virtual Compton process is obtained both in 
the form of a factorization theorem and an operator product expansion. In 
deeply virtual Compton scattering, they can be probed through distinct an- 
gular dependence of the cross section. 
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I. INTRODUCTION 



From the point of view of quantum field theory a nucleon is fully described only if 
one knows the matrix elements of all possible quark and gluon operators involving the 
nucleon state. Nevertheless, progress is possible provided one can obtain the matrix element 
of operators with a clear physical interpretation. Twist-two operators give the leading 
contribution in appropriate hard processes, are relatively simple, and are more accessible 
to experimental measurement. The states between which these operators are sandwiched 
may be of equal or unequal momenta; the former situation is familiar from well-investigated 
processes like deeply inelastic scattering. The latter has also been investigated over the years 
IjJ. However their importance has been understood only recently. For instance, knowing 
certain off-forward twist-two matrix elements allows for extraction of the quark and gluon 
orbital and spin contributions to the nucleon spin M. The class of parton distributions, 
known as off-forward (off-diagonal, non-forward) parton distributions defined from these 
off-forward matrix elements has generated considerable contemporary interest ^- [11 , 15 , 13 



In simple physical terms, a parton distribution, whether forward or off-forward, arises from 
removal of the parton from the nucleon by a hard probe and its subsequent return to form 
the nucleon ground state further along the light-cone, or a similar process. 

The class of twist-two operators which depends on parton helicity change shall be the 
subject of this paper. It is now well known flTl| , |TE| , |r7| that there is one forward chirally odd 
twist-two proton structure function, known as hi(x,Q 2 ), measurable in, for example, the 
Drell-Yan process. It is defined as the light-cone correlation of quark fields weighted by a^ u , 

^e iXx (PS'\ij(-l\n)a^^(lxn)\PS) = hx^UiPS'^UiPS) + ... . (1) 

Z7T 2 2 

In a helicity basis, wherein spins are measured along the particle's momentum and En is 
diagonal, S and S' differ by one unit of angular momentum. Since chirality and helicity co- 
incide for massless quarks, hi is a chiral-odd quantity. To give it a probability interpretation 
requires using a transversity basis wherein £j_7o is diagonal. This gives h\ the necessary 
probabilistic interpretation (for states of equal momenta): it is the probability to find a 
quark polarized along the transverse polarization of the nucleon minus the probability to 
find the quark polarized in the opposite direction. There is no gluonic helicity-flip distri- 
bution for obvious reason - a transverse gluon flipping its helicity leads to a change of two 
units of angular momentum and angular momentum conservation forbids this for a spin-| 
hadron. For hadrons with spin > 1 there is no such restriction and some years ago Jaffe 
and Manohar identified a leading twist gluonic structure function A(x,Q 2 ) which can 
be measured from a transversely polarized target like the deuteron. 

New Lorentz structures emerge if one allows for off-forward matrix elements, leading 
to generalizations of the above mentioned helicity changing structure functions. Recently, 
Collins et al. have suggested measuring the helicity-flip quark distributions in vector meson 
production [|l^,^]. The evolution equation for these distributions has been derived by 
Belitsky and Miiller ||11|| . Diehl et al. |21| have noticed that the distribution in angle between 
the lepton and hadron planes in deeply virtual Compton scattering (DVCS) contains valuable 
information about the helicity structure of the nucleon-photon amplitudes. They point out 
that photon helicity flip is possible even with a spin-^ target because gluons in the off- 
forward scattering can transfer two units of angular momentum. This, of course, requires 
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the existence of the gluon helicity-flip distributions in the nucleon. Indeed, for off-forward 
matrix elements one does not need a state of spin > 1 to accomodate gluon helicity flip. 

This paper is intended to present a comprehensive study of leading-twist helicity-flip 
off-forward distributions in the nucleon. A systematic counting suggests that there are four 
such distributions: two related to gluon helicity flip and the other two to quark helicity flip. 
In each case, a distribution can be defined depending on whether the nucleon's helicity is 
flipped or not. We derive the leading- logarithmic evolution of these distributions, although 
in the quark case, the result was already obtained by Belitsky and Miiller |15j . In the forward 
limit, the evolution of the quark distributions reduces to that of hx(x) as calculated by Artru 
and Mekhfi ||16|| . The evolution for the gluon distributions reduces to that of A(x), which 
has not appeared in the literature before. Note that the quark helicity-flip distributions do 
not mix with any gluon ones, and vice versa. This is quite significant because, for the first 
time, we have a parton distribution that can serve as a unique signature of the gluons inside 
the nucleon: gluonic effects cannot be mocked up by any kind of constituent quarks, and 
they cannot be evolved away by recklessly evolving down the momentum scale! 

We study measurement of the helicity-flip gluon distributions in general two-photon 
process. The photon helicity-flip Compton amplitude is calculated in terms of the gluon 
helicity-flip distributions. In the forward limit, we recover the result obtained in [T^|. How- 
ever, our result is more general. In the language of operator product expansion, we obtain 
the leading-order coefficient functions of a class of gluon operators with total derivatives. 



According to [pi] , the helicity-flip gluon distributions generate distinct angular distributions 
in the DVCS cross section. 

The presentation of the paper is as follows. In Section II, we enumerate the independent 
helicity amplitudes for the quark-nucleon and gluon-nucleon sub-processes. Subsequently 
new helicity changing distribution functions are motivated and defined. In Section III the 
leading logarithmic evolution of these functions is studied. Section IV contains a calculation 
of the Compton amplitude for photon helicity flip scattering. This amplitude vanishes at the 
tree level and requires at least one quark box (plus permutations) to be non-zero. Section 
V presents the DVCS cross section that depends on the gluon helicity-flip distributions. We 
conclude the paper in Sec. VI. 



II. HELICITY-FLIP PARTON DISTRIBUTIONS: 
COUNTING AND DEFINITIONS 

We shall, in this section, enumerate the complete set of off-forward quark and gluon 
distributions at the twist-two level. Helicity-flip ones will emerge through the counting 
and will be the focus of this paper. As usual, p M and n M are two light-like vectors with 
p 2 = n 2 = and p ■ n — 1. The momenta and spins of the initial and final nucleons are, 
respectively, P, S and P', S' . The momentum transfer A M = P ,fl — P M has both transverse 
and longitudinal components. It is convenient to define a special system of coordinates 
wherein P M = (P' + P) fl /2 is collinear and in the z direction, 

p> = p» + (M 2 /2)n >1 , 

A M = -2f (p" - (M 2 /2K) + , 

M 2 = M 2 - A 2 /4 . (2) 
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The initial nucleoli and parton have longitudinal momentum fractions 1 + £ and x + £, 
respectively. 

From dimensional reasoning, the leading order contribution to a given hard process must 
involve the minimum number of independent parton fields which, for QCD quantized on 
the light cone, is two (-0+ for fermions and A± for gluons). Therefore, one need consider 
only the matrix elements of bilinear operators at two different points on the light-cone. 
In the kinematic region x > £, one has the simple interpretation that the first operator 
extracts a certain type of parton from the nucleon and the second replaces it further along 
the light-cone. Let H, H' denote the respective helicities of the initial and final nucleon and 
h, h! the helicities of the parton extracted and replaced. The helicity amplitude Anh,H'h' 
must obey Awh',Hh — Anh,Wh' (time-reversal invariance), and A-h-h-h'-h 1 = Auh^ww 
(parity invariance). For a purely collinear process there is no preferred transverse direction; 
rotational invariance around the collinear axis requires helicity to be conserved, H + h' = 
H' + h. However, non-zero transverse momentum of the scattered nucleon or parton means 
that, while the total angular momentum will of course be conserved, helicity conservation 
will not necessarily hold. The difference is, of course, absorbed by the orbital motion of the 
scattered pair. 

For quarks it is readily seen that a set of independent amplitudes is provided by the 
following: An n, Ai_i Aii_i_i An Ai 1 _i 1, and A±_i _i 1. The familiar 

22'22 2 2'2 2 2 2' 2 2 2 2'2 2 2 2' 2 2 2 2' 22 

distributions fi(x,Q 2 ), gi(x,Q 2 ), and hi(x,Q 2 ) are linear combinations of the first three in 
the forward limit ]17|]. One unit of orbital angular momentum, made available by one power 
of Aj_, allows for three additional amplitudes. A complete set of off- forward leading- twist 
quark distributions is given below: 

I toi^p,^ 9 (-~An) 7 ^ A*n)\PS) =H q {x^)U{P' S')^U{PS) 

+E q (x,OU(P'S>)-^U(PS) + ..., 

J ^ c **<p^|&(-iAn)7^ 

+Eq (x,OU(P'S')-^ w U(PS) + ... , 
f ^^ {P i S >\$ q {-\\n)a^ q {\\n)\PS) = H Tq (x^)U(P> S')a^U(PS) 

J ATI L I 

+E Tq ( X ,OU(P'S')^^U(PS) + ... , (3) 

where \pu] means antisymmetrization of the two indices and the ellipses denote higher twist 
structures which are outside the scope of the present discussion. The dependence of each 
distribution upon t = A 2 and Q 2 is implicit. In each equation, the first term represents an 
amplitude that survives the forward limit and the second term an amplitude that decouples 
(but does not vanish) in the forward limit. The definitions of H q , E q , H q , and E q are from 
Ref . 0] . The quark helicity-flip distributions Ht q and Et q are new; and they can be selected 
from the third equation by taking \i — + and v =_L. The above definitions complete the 
identification of all twist-two quark distributions. 

A few additional comments about the definition in Eq. (|3D are in order. First, for brevity 
we have not explicitly shown the gauge link between the quark fields. This link is always 
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present, except in the light-like gauge A + = 0. Second, by using time-reversal symmetry, all 
the distributions are seen to be real. Third, from taking the complex conjugate of the above 
equations, it follows that all the distributions are even functions of £. Finally, we can add a 
time-ordering between the two fields without changing the content. This follows from, 

T4(0)^ + (An) = Vt(0)V+(An) - 0(An°){Vi(O)^ + (An)} . (4) 

The second term is just a constant because it is an anticommutator of the independent (or 
good) components of the Dirac field separated along the light cone. Obviously, the constant 
does not contribute to the matrix elements. More elaborate but essentially equivalent proofs 
can be found in the literature [p2,23fl . 

We now turn to the gluon distributions. Only transverse gluons need be considered 
here because longitudinal ones are either dependent or gauge degrees of freedom, which 
lead to either higher twist distributions or gauge links. An independent set of nucleon 
gluon amplitudes is: A^^, Aix_i_i -4.i lj i_ 1 , Aii_ii, and Ai_ lt _x x . The 

familiar distributions G(x, Q 2 ) and AG(x, Q 2 ) come from the forward limit of the first two 
amplitudes. There is no equivalent of h\ for gluons since it is impossible for a nucleon to 
spin-flip by two units. All off-forward twist-two gluon distributions are defined below: 

\ I ^^( P ' S ^ F{m (-\ n )^\\n)\PS) 

= H g (x^)U(P>S')P^U(PS) + E g {x,t)U{P'S') Hups) + ... , 

I J ^(P'S'\F^(-±n)iFj\±n)\PS) 

= H g (x, $U(P'S')P^hsU(PS) + E g (x, OUiP'S' f^^ UiPS) + ... , 

I f ^ e ^/p's'\F^ a (--n)F^\-n)\PS) 
x J 2n 2 2 

= H Tg (x,Z)U(P'S') m U(PS) 

+ E T9 (x,0U(P'S') — 7 - U(PS) + .... (5) 

Here in the first two equations, {^v) means symmetrization of the two indices and removal of 
the trace, and in the third equation [fia] and [p(3\ are antisymmetric pairs and (• ■ ■) signifies 
symmetrization of the two and removal of the trace. These operations are essential since the 
product of operators must transform as irreducible representations of the Lorentz group. 

The distributions H g , E g , H g , and E g have been introduced before Rll]. Their evolutions 
and mixing with quark distributions have also been worked out. They play an important 
roles in electro-meson production in the small x region P,p4| . The helicity-flip distributions 
Et 9 and Ht 9 are new. They can be selected from Eq. (||) by taking \i — v = + and 
a,j3 =_L. The angular momentum conservation requires presence of one unit of angular 
momentum (A^) when the nucleon helicity is flipped {Ht 9 ) and two units (A^A^) when it 
is not {Et 9 ). Hence both decouple from the matrix elements in the forward limit although 
the distributions themselves do not vanish. 
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In a sense, the gluon helicity-flip distributions are the "cleanest" among the class of gluon 
distributions since they are forbidden to mix with quark distributions by angular momentum 
conservation. This, in fact, was why Jaffe and Manohar [Tj| had proposed using the A(x, Q 2 ) 
distribution as a probe of "exotic gluons" in a spin J > 1 nucleus. Nuclear binding or pions 
would not contribute. However there, unlike here, there is a strong suppression on the 
magnitude of the distribution due to the small size of the nuclear interaction relative to a 
typical hadronic mass scale because "exotic gluons" can only be generated by the nuclear 
interaction (for an estimate of A(x, Q 2 ) see Ref.( ||25|| ) ). In the nucleon, the off-forward 
helicity-flip gluon distributions can be as large as other gluon distributions for a reasonable 
size of x. It is interesting to note that the n=2 moment of the gluon distributions above can 
be directly expressed in terms of the matrix elements of the colour electric and magnetic 
fields in the nucleon. 



III. EVOLUTION OF HELICITY-FLIP DISTRIBUTIONS 



A parton distribution is necessarily defined at a given distance or momentum scale be- 
cause of ultraviolet divergences. Physically, the scale can be related to the kinematic vari- 
ables of the particular hard process under consideration. In this section we study, using 
momentum space Feynman diagrams, the leading-logarithmic evolution of the helicity-flip 
parton distributions defined in the previous section. While the method is straightforward 
in principle, there are important subtleties related to the gauge dependence of the calcu- 
lation and end-point singularities. In principle any gauge choice is permissible but the 
light-cone gauge is naturally preferred for several reasons: the fewest number of diagrams 
need be calculated, path-ordered exponentials are absent, the light-cone dynamics of partons 
is transparent, and ghosts are absent. However, there is a price to be paid because the light- 
cone gluon propagator, when imbedded in a loop, leads to singularities in the end-points 
of integrals whose interpretation is ambiguous. This is sufficient reason to do the (longer) 
calculation in a covariant gauge. Therefore, in the following, we shall use the Feynman gauge 
and treat separately the quark and gluon helicity-flip distributions. At leading-logarithmic 
order, any ultraviolet regulator is as good as any other; so we impose an ultraviolet cut-off 
on the momentum integrations. 



A. Evolution of the Quark Distributions Et and Ht 



The result in this section has been obtained before in Ref. pT| . Here for completeness 
we present our calculation in a different form. 

From the definitions of the quark distributions Htq and Et q in Eq.(|3]), it is clear that 
the two will evolve identically with the leading component of the operator ipa^ip. Selecting 
only the leading twist part, it is therefore convenient to define, 

Hx^Q 2 ) = n,e u J ^^P'S'm-^Xn^e-* 9 ^ (6) 

where e M is a unit vector in a transverse direction. The path-ordered integral has been 
reinstated in the above. The renormalization scale Q 2 is the cut-off for the momentum 
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components of the fields. J-(x, £, Q 2 ) can be diagrammatically represented as in Fig. la. Now 
imagine a slight increase of Q 2 to Q 2 + 5Q 2 , revealing a deeper level of hadronic substructure. 
Additional diagrams contributing to J-{x, £, Q 2 + SQ 2 ) are shown in Figs, lb-le. 
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FIGURES 



(a) 




(b) 



(c) 



(d) 



(e) 



FIG. 1. Feynman diagrams for evolution of quark helicity-fiip parton distributions. 

Feynman-like rules for these diagrams can be derived in a rather straightforward way: 
the product of bare Heisenberg fields in Eq.(|6]) can be brought under the time-ordering 
symbol and a perturbation expansion follows from expanding out the exponential containing 
the QCD interaction terms exp(i / £i nt ) as well as expanding the path-ordered exponential. 
Subsequently a collinear expansion is made: the lines entering or leaving the hadron blob 
are always collinear, reflecting the fact that the internal relative momenta of partons in the 
hadron are much less than the momentum of the hard probe. The transverse momenta of 
the other lines is bounded by Q 2 . 

The Feynman expressions for fig.lc and Id are, 



5F c {x,£,Q 2 



SF d (x,t,Q 2 



(q — k) ■ n + ie 



iD aP (q - k) 



xTr 
1 



-a pu {-igt a n a )a p \-igt a ^)iS F {k + -A) 
8 2 



2 n ^J dV J (27^ 



d 4 k 



8(x - y) 



(k — q) ■ n + ie 



iD a/3 (q - k) 



(7) 



xTr 



-o pu {-igt a n a )a p \-igt a ^)iS F {k + -A) 
8 2 



p [p e x] F(y,£,Q 2 ) . (8) 



In the first of the above two equations, q p = yp p and in the second, q p = xp^ 1 . Two other 
diagrams are obtained by reflection and yield identical expressions under £ — > — £. Finally, 
the quark self-energy diagram in fig.le (and its reflection), together with the wavefunction 
renormalization, yield the combined expression below after calculation, 



D Q T{x,^Q 2 ) a s {Q 2 



DlnQ 2 



2tt 
-0{& 



-Cr 



X 



6(x-0 
x dyx - £ 

-i y y-Z 



1 dy x — £ 



V v-t 

- 0{-t - 



+ 0(x + £) 



dy x + £ 

y y + £ 



X 



dyx + i 



-i y y + i 



Hy,ZiQ 2 ) 

y — x + ie 



(9) 



where, 
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d 



a s (Q 2 

D\nQ 2 d\nQ 2 2vr 



3 + P dy 
2 J% y — x — ie 



+ 



dy 



e, y — x — it 



(10) 



The above result agrees with that obtained by Belitsky and Miiller. 
It is useful to take moments of T{x, £, Q 2 ). Define, 



F n (.€, Q 2 ) = f_ dx x n ~ l F{x, i, Q 2 ) (n > 2) 



(11) 



where n =even (odd) moments are charge conjugation even (odd). Then, a calculation leads 
to the following evolution equation for the moments, 



dT n 
rflnQ 2 



a s {Q 2 
2tt 



-C, 



1 2 J 1 1. 



where, 



--ZS{n))F n + 2 E (---)(20 2i ^n- 

Z J t=l,2,.. Zi 72 



n i 

£1 * 



2/ 



(12) 



(13) 



We shall now interpret the evolution equation in terms of operator mixing. To this end, 
define, 



_ <-> A^i+l <-» Mn— 1 

= id" 1 ■ ■ ■ i&^ij) iD ■■■iD a^Tp . 



(14) 



where all \x indices are symmetrized and iD= [iD — iD]/2. Using translational invariance, 
it is easy to see that, 



n 



in 



n—2i 



(15) 



As one can see from Eq . (|12|) , the operators O^i^ belonging to same n but different i mix. 
Taking an appropriate linear combination of these, or equivalently, diagonalizing the mixing 
matrix, is not difficult. One may establish recursion relations using Eq.(^) to finally arrive 
at, 



ip, = (i _ 2S(n) , o. , 

dlnQ 2 2tt V2 K n 



r n — 1 ] 

\. 2 J ( ]\% c 2n — 2i—l 



E 

i=0 



T(n-i + ± 



(n-2i-l)! i\ T(f) 



2> (T)Hi---IMiV 



(16) 
(17) 



The coefficients in Eq.([l7]) are those of the Gegenbauer polynomials, C%-i(x)- This can 
be traced to the fact that these polynomials are essentially the Clebsch-Gordan coefficients 
which occur in the light-cone expansion of operator products that transform irreducibly 
under the conformal group K5 . 
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B. Evolution of the Gluon Distributions and H? 



The evolution of the gluon helicity-flip distributions can be studied in the same way as 
for the quark. For convenience, we define, 



F(x,£,Q 2 ) = n^e^e'^ 

dX 
— < 

2tt 



x / _ e ^(P / 5 , |F^(-iAn)e" l9 fi /2dan - A(an) F^(iAn)|P5). 



(18) 



where e a and are two unit vectors in the transverse directions. We are interested in the 
change of F under the change of the momentum cut-off Q 2 . Some representative Feynman 
diagrams are shown in Fig. 2. 






(a) (b) ' (c) 

FIG. 2. Representative Feynman diagrams for evolution of gluon helicity-flip parton distribu- 



tions. 



As we have mentioned before, the gluon helicity-flip distributions do not mix with any 
quark distributions. For x > £, the evolution equation reads, 



D Q F(x,£,Q 2 ) = a s (Q 2 ) /•! dy p fx £ e 
DlnQ 2 2vr J x y \y'y'y. 



F(2/,£,Q 2 ) 



(19) 



where, 



D, 



d a s (Q 2 



DlnQ 2 dlnQ 2 2tt 



-c A 



11 n 



f 



+ 



dy 



+ 



6 3Ca h y — x — ie J-£ y — x — ie 



dy 



(20) 



and, 



' x (l_^2)2 ^ x 2 _£2 ' l_ x + i£ 

For — £ < x < £, the equation is 
D Q F(x,£,Q 2 ) = a s (Q 2 ) ( p dy pl fx £ _ f* dy p , (x £ e 

where, 



(21) 



2vr U y \y y y, 



. . , ,F(j/,e,g 2 ) , (22) 

i y \y y y 



(x 2 - 



2? 



+ 



(x- 0(1 + 1-x + ie 



(23) 



10 



And finally for x < — £, the equation is the same as that for x > £, except J^ 1 — > — f^ v 
In the forward limit, the evolution equation reduces to that for A(x), 



dA(x) 
dlnQ 2 

where the evolution kernel is, 



a s (Q 7 
2tt 



dy 

x y 



P 



x 



My) 



P(x) 



2G, 



X 



1-x) 



+ 



( HC A _ nj_ 
V 6 3 



6{x-l) . 



(24) 



(25) 



Here the + prescription is standard |27|]. The above result, as far as we know, is new. 

It is instructive to look at the evolution in operator form. Define the n = even moments 
(the odd moments vanish because F(x) is antisymmetric in x as one can easily check from 
the definition), 



The evolution equation becomes, 



(n > 2) 



(26) 



dF r 



dlnQ 2 



a s (Q 2 
2vr 



-Ca 



11 



n 



f 



6 2C, 



2S(n) F n + £ 



i=i 



4z - 2 4z + 2 1 
+ ■ 

z 



+ - ) (20 2i F„_ 2i 



Define a tower of twist-two gluon operators, 



O 



n,2i 

Then it is easy to see 



n 



/'I 



n,2i 



|P5) = (20 2l F n _ 2 



(27) 



(28) 



(29) 



Hence the mixing of the different moments of the gluon distributions reflects the mixing of 
the twist-two operators of same spin and dimension. Define a new basis of operators in term 

5 

of the Gegenbauer polynomials C^_ 2 (x) combinations, 



r re— 3 1 



E 



r(n — z + 



4=0 (n-2z-2)!z!r(|) 
Then the evolution of O n becomes diagonal, 



2 1 ssyn—HnCcP 



(30) 



a s (Q 2 )„ (11 



dlnQ 2 



2tt 



6 3C F 



2S , (n) O 



(31) 



Again, the above simplification is due to conformal symmetry. However, beyond the leading- 
logarithmic order, the conformal symmetry is anomalously broken by quantum corrections 
2|. 
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IV. PHOTON HELICITY-FLIP COMPTON AMPLITUDE 



The gluon helicity-flip distributions, H Tg (x,^) and E Tg (x,£), are basic properties of 
the nucleon, at par with the other distributions, namely, H g (x,(,), E g (x,£,), H g (x,£,), and 
E g (x, £). Because of angular momentum conservation, the gluon double- helicity flip distribu- 
tions do not mix with quark distributions, making their isolation and possible measurement 
relatively cleaner. One would like to know which hard processes probe Hj< g and -Erg- The 
general Compton scattering involving two photons offers one possibility, perhaps the sim- 
plest. Diffractive vector meson production from a deeply- virtual photon is another |19| , |24 




(a) 



V 





(b) (c) (d) 

FIG. 3. Diagrams for photon helicity- flip Compton scattering. 

Fig. 3a illustrates the general two photon process. A highly virtual photon with mo- 
mentum q + A/2 and transverse polarization is incident upon the nucleon. A second photon 
of momentum q — A/2 is detected, and the recoiling nucleon emerges intact. The scattering 
amplitude is, 



= i J d*z e- iq - z (P'S'\T 



r{-\z)n\z) 



\PS). 



(32) 



A convenient set of kinematic variables is obtained by choosing and = |(P + P')^ 1 
to be collinear and in the 3 direction. In terms of n M , we expand the "average" photon 
momentum, 



<f = - XbV » + (Q 2 /2x B )n» . 



(33) 



The transverse components of A M are assumed to be of the order of the nucleon mass, i.e. 
much smaller than the hard momentum Q. In the collinear approximation, the momenta 
of the incoming and outgoing photons are, respectively, q — £p and q + £p. The initial and 
final state nucleon's momenta are (1 + £)p M and (1 — £)p M , respectively. Symmetrizing the 
diagrams in this manner makes it possible to get simpler expressions by exploiting crossing 
symmetry. 
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There is no tree-level coupling between the scattering photon and the gluons in the target. 
However, there is a contribution at 0(a s ) from the box-diagram and its permutations, shown 
in Figs. 3b-d. The momenta of the gluons going outward and returning to the nucleon blob 
are, respectively, (x+Op^ and (x—^p^ with — 1 < x < 1 . The loop momentum is expanded 

as, 

k» = (k- n)p^ + (k ■ p)n» + k^ . (34) 

Rather than use Feynman parameters to solve the integral with four products in the denom- 
inator, we systematically exploit crossing symmetry (£ — > — £) to reduce the integral to a 
much more manageable form with three product terms. Dimensional regularization is used 
for the divergent graphs. The infrared singularity cancels in the sum of the three diagrams, 



as it must. Many details of the calculation are similar to those in Ref. [£S 
Our final result for the photon double helicity-flip amplitude is, 



where, 




1 + —5 i. l n —5 !L 



n a n p T^ , (35) 



T v»ap = i I ™ e i^(p's'\F^(-\)ir®(-n)\PS) . (36) 
J ^"Tc 2 2 

In the above, q sums over the quarks circulating in the loop and e q are their electric charges. 

Eq. ( PoD is the convolution of a perturbatively calculable part and a soft part reflecting the 

nucleon's composition, and is in the form of a factorization theorem at the lowest order 

of QCD perturbation theory. To completely justify the use of the soft part, one needs 

to consider all contributing Feynman diagrams with the hard quark loop connecting with 

arbitrary number of gluons to the nucleon blob. Any choice of gauge may be made of course, 

but the final result will be gauge independent. It is easiest to work in the light-cone gauge 

and with simply two physical gluon fields. In covariant gauge it would be necessary to show 

that these additional gluons are summed up in the path-ordered integral. 

As a check on our calculation, we may consider forward scattering £ = on a target with 

spin J > 1. In the limit of xb — > oo, the perturbative part in Eq.(^) can be expanded: 

T ,u = _^ (y e 2) y — — f 1 dx x n - l T^n a n p 




n=2 n + 2 (Q^ 

x {PS\F^ 1 iD ■■■iD F U ^\PS) (37) 

The above coincides with the result in Ref. |18j except for an overall sign. In the general 
case, we can convert Eq . (|3~5l) result into a generalized operator production expansion with 
derivative operators, 

i J d A ze^ z J v {- Z -)J^) = ^ ^ ^ 2 "' € 

x £ {n-2i)O n *.^ v + ... (38) 

i=0 
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where Offi.n n uu * s defined in Eg. fl2"8|) . We would like to emphasize again that only the \iv 
symmetric and traceless terms are included in the above equation. 

V. DVCS CROSS SECTION WITH GLUON HELICITY-FLIP DISTRIBUTIONS 

The Compton amplitude in the last section can be used to obtain the cross section for 
deeply virtual Compton scattering. In DVCS, the final photon is real and hence one has the 
constaint xb = £■ In this section, we calculate the cross section in this special kinematic 
limit. 

We choose the kinematic variables as those used in Ref. ||: k = (u,k) and k! = (u/, k') 
for the four-momenta of the intial and final electrons, P = (M, 0) and P' = (£", P') for the 
initial and final momenta of the nucleon, and q' = (u', q') the momentum of the final photon. 
The differential cross section is, 

da = \T\ 2 dT , (39) 

where T is the invariant T matrix of the scattering and dT is the invariant phase space factor. 
Depending on choice of independent kinematic variables to characterize the differential cross 
section, dT takes different form. If one uses the scattered electron's energy and solid angle, 
and the scattered nucleon's solid angle, 

1 P' 2 

dT = — — -——uj'duj'dVLe'dVLp,-— — , (40) 

32(2tt) 5 cjM y \P'(u + M) - qE 'cos0| ' v ' 

where is the angle between q and P', and the sum over two possible solutions of \P'\ is 
implicit. On the other hand, one can also use the standard Q 2 and xb (or s), t = A 2 , the 
t-channel momentum transfer, and the angle between lepton and hadron planes [pi| , pO|| . 

The invariant T matrix consists of two parts. The first part comes from the Compton 
scattering, 

r x = -juW-fuW^T^fr , (4i) 

where u, u are the spinors of the lepton and e is the polarization of the emitting photon. 
The Compton amplitude T^ v contains both photon helicity-flip and non-flip contributions, 

rp ry-iAA=0 _i_ T"iAA=2 (AO\ 

where the first term is given by Eq.(4) in Ref. and the second term is from Eq. (|35|) . 



47r \Y / \x — £ + ie x + £ — ie 



P^iA al a u ^ 
H Tg (x,0U(P'S') ^ U(PS) 



P^Ai 4 V A P ? 
+ E Tg (x,0U(P'S>) u ]f IH'S) 
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Air 

x 



x - £ -He x + ^-ieJ AM 
[H Tg (x,OU(P')(A^ + AM - A±) iU{P) 



+^E Tg (x,OU(P') + - 9? A±)+ ^(A^AI - i<TAi) ) C/(P) 

The second part of the T matrix comes from the Bethe-Heitler process, 



(43) 



% = -e 3 u(k') 



1 



A — m e + ie 



7 M + 7%7 



1 



+ A — m e + ie 



u{k)±- 2 (P>\JM\P) , (44) 



where m e is the mass of electron and will be ignored for the following discussion. The elastic 
nucleon matrix element is, 



(p'\jM\p) = u{p') 



7At F 1 (A 2 )+F 2 (A 2 )^ 1 



U(P) , 



(45) 



where U, U are the nucleon spinors and Fi and F 2 are the usual Dirac and Pauli form factors 
of the nucleon. 

We are interested in only the leading contribution to the cross section from the helicity- 
flip gluon distributions. This comes from the interferences between the helicity-flip and 
non-flip Compton amplitudes and between the former and the Bethe-Heitler amplitude. 
The first interference yields, 



e° 



(46) 



//j-AA=0^* q-AX=2 _|_ q-AX=0 ^j-AX=2\* _ 

where the lepton tensor £y U c can be found in Ref. ||. The hadron tensor is 
^c = ^(AlAl-i,rAi)^(E^) 

x J2 e 2 q Re f 1 dxa{x) dx'a*{x') \H q (x, Z,)E Tg {x' , f ) - H Tg (x', £)E q {x, £)1 • (47) 



The tensor structure (A^A^ — ^g^ u A 2 ± ) signals a cos 20 term in the cross section, as was 
noted in [21|. There is, of course, also a gluon helicity non-flip term but it enters at one 
higher power of a s . 

The interference between the double-helicity-flip Compton and Bethe-Heitler amplitudes 

is. 



7i AA=2 ~r 2 + r^ A=z r 2 



AX=2n~* 



2/02 ^^(/i^a ) 



A 2 g 



(48) 



where £^ u > a depends on electron kinematic variables and can be found in [Ij. The nucleon 
structure dependent part is, 
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jj{hv)c 



A'lAl - -ft A 2 




dxa(x) 



x 



(Fx + F 2 



' A 2 
H Tg {x,£) + —Eh 



+ (A^r + a^j 




+ (iqi? Ts (x,£) - F 2 H Tg {x,i)) — 



dxa(x) 



(49) 



Here the presence of A in F m and A^A V in H^ a can give rise to a distinct cos 30 terms 
in the cross section [pT|. To obtain the latter, one just multiplies the lepton and hadron 



tensors together to get the square of the T matrix. Because of its length, we omit the final 
expression. 

A more direct way to see the angular dependence of the cross section is to use the 
formulas derived in the center-of-mass frame in Refs. According to these works, all 

A A' 

one needs is the hadron helicity amplitude M H ' H ,, where A, A and H, H are the initial and 
final photon and nucleon helicities, respectively. The helicity-flip nucleon amplitude M^ 1 ^, 
clearly is just the helicity-flip Compton amplitude, 



M, 



-1,1 



h,h' 



e a (-l)e*J+l)Tf X: 



=2 • 



(50) 



This can readily be evaluated using Eq.(f|3|) by substituting in the appropriate Dirac spinor 
for the nucleon helicity states. According to [pi], certain angular weighted cross sections 
can be used to make a direct extraction of the above amplitude. 



VI. SUMMARY AND COMMENTS 

In this paper, we have presented a number of new results related to the helicity-flip 
off-forward parton distributions. First, we enumerated systematically all leading- twist off- 
forward parton distributions for a nucleon: six for the quark parton and another six for 
the gluon. Four of these distributions, two each for the quark and gluon partons, involve 
parton helicity-flip. Second, we derived the leading-logarithmic evolution equations for these 
helicity-flip disitributions. In the forward limit, our result agrees with the known kernel 
for hi(x) while the kernel for A(x) is new. Third, we obtained the photon helicity-flip 
Compton amplitude in terms of a tower of gluon operators with total derivatives. Our 
result may be obtained from the known forward case by using the conformal symmetry of 
QCD. Finally, we compute the leading DVCS cross section which depends on the gluon 
helicity-flip distributions. 

We have emphasized the unique role played by helicity-flip distributions in characterizing 
the properties of the nucleon. If one askes for a clear experimental signal of existence of 
gluons in the nucleon, the helicity-flip Compton amplitude would serve the purpose. Without 
the vector gluons, it would be at least power suppressed in the high-energy limit. Of course, 
the helicity-flip gluon distributions can also be measured in vector meson production ||TU| . 



The size of the helicity-flip distributions should be similar to the usual helicity-dependent 
parton distributions - there is no extra suppression in the soft physics to curb helicity flip. 
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